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1. Introduction 

Of the 123 planets found so far around 108 main-sequence stars, there are 13 multiple 
planet systems. The fourth extrasolar planet discovered was the 14.65 day period planet 
around Cancri (= rho Cancri A = 55 Cancri 1= HD 75732 = HIP 43587 = HR3522) 
(Butler et al. 1997). More recently, the orbits of an additional outer planet and possible 
intermediate planet have been reported (Marcy et al. 2002). In this Letter, we announce 
the discovery of a fourth inner planet in the Cancri system with radial velocity semi- 
amplitude (K) of 6.7 m s~^ , an orbital period of 2.808 days, Msin i of 0.04 Jupiter masses, 
and a probable actual mass of 0.056 Jupiter masses (17.7 Earth masses). 

Cancri appears to be a super metal-rich (CSV, V = 5.95, [Fe/H] = -1-0.27), but 
otherwise normal main sequence star with a V=13 visual binary companion. There is no 
detectable circumstellar disk around Cancri (Jayawardhana et al. 2002; Schneider 2001). 
The metallicity is significantly higher than the mean of nearby field stars, but is in accord 
with the trend of close-in giant planets being found preferentially around metal-rich solar- 
type and cooler stars (Taylor 1996; Gonzalez & Vanture 1998; Feltzing & Gonzalez 2001; 2; 
Laws et al. 2003; Heiter & Luck 2003). Mass estimates for p^ Cancri range from 0.88 - 1.08 
Mq (Laws et al. 2003; Fuhrmann et al. 1998; 3; Santos et al. 2003). We adopt a mass of 
0.95 Mq ± 0.08 (AUende-Prieto 2004, personal communication). 

2. Observations and Results 

2.1. HST Astrometry 

This Letter made use of reanalyzed public domain data (McArthur et al. 2005, in prepara- 
tion) obtained from the Hubble Space Telescope {HST) (McCrath et al 2002; McCrath et al 2003) 
to provide an estimate of the inclination of the long period planet, p^ Cancri d. For the pa- 
rameters critical in determining the mass of p^ Cancri d, we find a parallax, iTabs = 79.78 ± 
0.3 mas, a perturbation size, a= 1.94 ± 0.4 mas, and an inclination, i = 53° ± 6.8° This 
inclination was derived from a small arc of orbit coverage in the limited HST data set. In 
any case, an inclination of less than 20° can not be considered because of the very large 
perturbation that could not have been missed, even by the small existing sample of HST 
astrometric data. 
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2.2. Radial Velocity 

We used the Hobby-Eberly Telescope (HET) to make many, high precision radial veloc- 
ity observations of Cancri over a short time period in preparation for anticipated additional 
HST astrometry, knowing that these observations would also be able to uncover any addi- 
tional planetary companions with radial velocities greater than 4 m and orbital periods 
of a year or less. The HET High Resolution Spectrograph (HRS), used in conjunction with 
an I2 gas absorption cell, can give routine radial velocity precision of 3 m s~^ for high signal- 
to- noise observations (Cochran et al. 2004). By taking advantage of the queue-scheduled 
operation of the HET, we recorded well over 100 radial velocity measurements in a span 
of about 190 days, shown in Table 1. This "snapshot" of the Cancri system promised 
the best definition of the orbital parameters of the two known short-period inner planets, 
because it would be relatively free of time- varying changes, not only in the planetary system 
but also in the facilities used to observe the system. 

The long period planet (p^ Cancri d) could clearly not be modeled by the HET data 
alone, but because of its velocity amplitude it was a very critical part of the initial modeling 
of the triple-Keplerian orbit. Three sets of radial velocity data from Lick (Marcy et al. 2002), 
ELODIE (6) and HET were fit with a triple-Keplerian orbit using Gaussfit (5), employing 
the techniques to combine RV from different sources (Hatzes et al. 2000) producing the orbit 
of the outer planet shown in Figure 1. This orbital fit remains problematic because of some 
discontinuity between the data sets for Cancri c (P=44d), and the sparse early Lick data. 
It is clear that the extrema of the HET data are not fit as well as the central portions. The 
orbit derived for this long period planet was then used as a constant in further modeling. 

p^ Cancri c was discussed (Marcy et al. 2002) as a signal that could possibly be caused 
by rotation modulation of inhomogeneous surface features. The stellar rotation period 
of 36-42 days was uncovered from periodicities in the Ca II emission (Henry et al. 2000; 
Noyes et al. 1984). However, the planetary hypothesis is more likely because Cancri is a 
very inactive star in which stellar activity induced RV efi'ects generally do not exceed 3.0 
m s~^ . The amplitude (Marcy et al. 2002) of the signal for this intermediate planet was 
over 12.0 m s~^ . The phased HET data showing p^ Cancri c are shown in Figure 2. Figure 3 
shows p^ Cancri b, the 14.66 day period planet, the original planet found around p^ Cancri . 

Analysis of the residuals of the HET data, after removal of the best overall fit to the 
orbits of the three known planets to all of the radial velocity and astrometric data, revealed 
a velocity perturbation with a period of 2.808 days. Simultaneous modeling of p^ Cancri b, 
p^ Cancri c, and p^ Cancri e with p^ Cancri d as a constant, showed an object in a low 
eccentricity orbit with a period of 2.808 ± 0.002 days. The radial velocity semi-amplitude 
of this orbit is only 6.67 ± 0.81 m s~^ , which is above the expected internal radial velocity 
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" jitter" variability for this star. The 2.808 day radial velocity signal is coherent over the entire 
time-span of the HET data, which would not be expected for stochastically excited stellar 
surface sources. The HET residuals, phased to the 2.808 day period with the orbital solution 
over plotted arc show in Figure 4. In retrospect, the clear signal of this 2.808 day orbital 
period planet can be seen (along with the signal of the 44 day periodicity (p^ Cancri c)) in 
Figure 9 of Marcy et al. 2002, which shows the periodogram of the residuals of the Lick data, 
after the Cancri b and d have been removed. This gives independent confirmation that 
this radial velocity periodicity is intrinsic to the Cancri system, and is not an artifact of 
the observation or data analysis techniques used on either the HET or Lick data. The 2.808 
day periodic signal has a false-alarm probabihty of 1.7 x 10~^ in the HET data (as seen 
in FigTuc 5) and 1.1 x 10~^ in the Lick data, giving further strong evidence for its reality. 
We interpret this 2.808 day radial velocity periodicity as the stellar reflex motion due to an 
extremely low mass planet in orbit around the star. Properties of the combined four-planet 
orbital solution to all of the radial velocity data are given in Table 2. Mass parameters and 
limits are shown in Table 3. The HET data in time overplotted with the quad-Keplerian fit 
are shown in Figure 6. The mass function from the radial velocity solution, combined with 
our assumed stellar mass of 0.95 ± 0.08 solar masses, gives an Msin i of this innermost 2.808 
day planet of 0.045 ± 0.01 Mjup . If we assume that all of the planets in the p^ Cancri system 
are coplanar and adopt the 53° inclination of the outermost planet, then we compute a true 
mass for this innermost planet of about 17.7 Earth masses, close to the mass of Neptune. (If 
an inclination of 20° is assumed, an upper limit of 2.37 MneptuneIS derived.) This fourth 
planet, Cancri b, is the lowest mass extrasolar planet yet found around a solar-type star. 
We note that the lowest mass, and indeed the very first, extrasolar planets ever found are 
those around the millisecond pulsar PSR B 1247+1 22 l(Wolszczan & Frail 1992). 



3. Discussion and Conclusion 

The inclination of the Cancri e, assuming coplanarity of the system , indicates that 
this object is indeed a Neptune-mass object, and not a more massive planet viewed nearly 
face-on to its orbit. The presence of this very low-mass planet at such a small semi-major 
axis in a system of at least three other gas-giant planets provides an important testbed for 
models of planetary system formation and evolution. 

A crucial question is whether p^ Cancri c originally formed at something near the current 
observed mass, or whether it formed as an intermediate or much more massive gas giant that 
has lost a significant amount of mass during the evolution of the system. 

It is reasonable to assume that at least the three inner planets of the p^ Cancri system 
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were formed at roughly the same time, and that their dynamical evolution was closely linked. 
The most massive of these three planets, p-*^ Cancri b (P=14.66d), would have been the 
first to achieve a mass sufficient to open a gap in the disk and begin Type II migration 
(Ward 1996). If this happened before the inner Cancri e reached the critical mass for 
rapid accumulation of a gas envelope, then it could have been carried inward as part of 
the overall migration of the disk material without ever growing past the rocky core phase 
(Lin & Papaloizou 1986). Whether this inward tidal migration driven by Cancri b would 
sweep such a Neptune-mass planet into the stellar photosphere, or would leave it in its present 
short period orbit depends on the details of the end-process of planetary tidal evolution. If the 
planet formed at a mass near its present mass by the core-accretion model of planetary system 
formation (Pollack et al. 1996), it would be a rocky body which could have a substantial 
dense atmosphere. 

An alternative hypothesis is that Cancri e formed as a much more massive gas 
giant planet before it began its inward type II migration to its present orbital radius 
(Ida & Lin 2004). With the planet then within 0.05 AU of the parent star, it could have 
suffered significant heating from tidal interactions with the star (1). This tidal heating 
could have caused the radius of the planet to infiate, possibly overfiowing the Roche radius 
of the star and resulting in significant planetary mass loss (4; Baraffe et al. 2004). It is 
possible that the mass loss has stripped this planet of almost all of its original H2-IIe gas 
envelope, leaving just the rocky planetary "core". This loss of a gas envelope for very short 
period planets is supported by the observation of an extended H I, O I and C II atmo- 
sphere around HD209458b (Vidal-Madjar et al. 2003; 2004), which may be in hydrodynamic 
escape from the planet (Lecavelicr ct al. 2004). However, the discovery of transiting Jupiter- 
mass planets (OGLE-TR-113 and OGLE-TR-132) with orbital periods of less than 2 days 
(Bouchy ct al. 2004) may cast doubt on this scenario. Both of the OGLE planets are about 
one Jupiter mass and are only slightly greater than one Jupiter radius. This indicates that 
planets which are much more massive than p^ Cancri e can indeed survive at even shorter 
semi-major axes. 

We have shown that an intense campaign of directed observations of a planetary system, 
aided by previous observations from other telescopes, can uncover Neptune mass objects. 
While it will be years before astrometric missions (SIM and GAIA) will be able to uncover 
the astrometric perturbations orbits of all these planets, we have been able to use HiSTdata 
to place a preliminary limit on the inclination of the outer planet. 

The Hobby-Eberly Telescope (HET) is a joint project of the University of Texas at 
Austin, the Pennsylvania State University, Stanford University, Ludwig Maximillians Uni- 
versitat Miinchen, and Georg August Universitat Gottingen. The HET is named in honor of 
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its principal benefactors, William P. Hobby and Robert E. Eberly. These results are partially 
based on data obtained with the UCO/Lick Observatory Telescope, the ELODIE spectro- 
graph at Observatoire de Haute-Provence (OHP, France), the Apache Point Observatory 
3.5m telescope which is owned and operated by the Astrophysical Research Consortium), 
and the NASA/ESA Hubble Space Telecope (obtained at the the Space Telescope Science 
Institute, which is operated by the Association of Universities for Research in Astronomy, 
Inc). This paper is based upon work supported by the NASA under grants HST GO-09969, 
NAG5-13206 and NNG04G141G, and NSF under grant AST-9808980. We are grateful to 
Ed Nelan for providing the HST d&ta, and Mehssa McGrath, for supporting our further in- 
vestigation, and to W. Spiesman for his careful reading of this manuscript. We are grateful 
to our anonymous referee for their thoughtful comments. 

Facilities: McD Obs(HET), HST(FGS), UCO/Lick), OHP(ELODIE), ARC(Apache 
Point). 
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Table 1. Measured Velocities for Cancri 



Julian Date 


Radial Velocity 


Uncertainty 




m s~^ 


m s~^ 


2452927.984100 


48.3959 


4.42 


2452927.986693 


45.2119 


4.61 


2452928.973951 


32.9745 


4.33 


2452928.976555 


33.8595 


4.00 


2452929.976107 


10.1059 


4.01 



Note. — [The complete version of this table is in 
the electronic edition of the Journal. The printed 
edition contains only a sample.] 



Table 2. Quad-Keplerian Orbital Elements of Cancri 



Element 


Cancri e 


Cancri b 


Cancri c 


Cancri d 


( iTrn'l'fi 1 Mon f\c\ i H cnrQ 1 


9 ono _)_ n nno 


1 4 «7 _)_ n ni 




4c;i 7 4 _(_ 77 Q 


Epoch of Periastron T ^ 


3295.31 ± 0.32 


3021.08 ± 0.01 


3028.63 ± 0.25 


2837.69 ± 68.87 


Eccentricity e 


0.174 ± 0.127 


0.0197 ± 0.012 


0.44 ± 0.08 


0.327 ± 0.28 




261.65 ± 41.14 


131.49 ± 33.27 


244.39 ± 10.65 


234.73 ±6.74 


Velocity amplitude K (m ) 


6.665 ± 0.81 


67.365 ± 0.82 


12.946 ± 0.86 


49.786 ± 1.53 


Lick (m ) 


21.166 ± 1.31 








ELODIE (m s-^ ) 


2727.448 ± 2.42 








^0 HET (m ) 


10.745 ± 0.59 









^Add 2450000.0 to T 



Table 3. Cancri - Ivbiss Limits and Parameters 



Parameter 


Cancri e 


Cancri b 


p^ Cancri c 


p^ Cancri d 


a (AU) 


0.038 ± 0.001 


0.115 ± 0.003 


0.240 ± 0.008 


5.257 ± 0.208 


A sin i (AU) 


1.694e-6 ± 0.19^6 


9.080e-5 ± 0.12^5 


4 .695^5 ± 0.14^5 


0.195e-l ± 0.007^1 


Mass Fraction (M© ) 


8.225e-14 ± 2.33e-14 


4.64e-10 ± 0.17e-10 


7.151e-12 ± 0.54e-12 


4.874e-08 ± 0.38e-8 


Msinz (Mjup Y 


0.045 ± 0.01 


0.784 ± 0.09 


0.217 ± 0.04 


3.912 ± 0.52 


Msinz (Mneptune)^ 


0.824 ± 0.17 








Msinz (Mearth)'' 


14.210 ± 2.95 








M (Mjup f''' 


0.056 ± 0.017 


0.982 ± 0.19 


0.272 ± 0.07 


4.9 ± 1.1 


M (Mjup f^" 


0.053 ± 0.020 


0.982 ± 0.26 


0.244 ± 0.07 


4.64 ± 1.3 


M (Mneptune)'^''' 


1.031 ± 0.34 








M (Mearth)'^''^ 


17.770 ± 5.57 









^derived from radial velocity alone 

''derived from radial velocity and astrometry, using Msini/sini 

■^derived from radial velocity and astrometry, using m2^/(ml + m2)^ = a^/P^ 

^assumes coplanarity of the planetary system 
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Fig. 1. — The RV orbit of the long period planet (p^ Cancri d) is shown with the three data 
sets, Cancri b and Cancri c are subtracted from the observations. The rms of the HET 
residuals in this fit are 9.2 m s~^ . 
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Fig. 2. — The RV orbit of the 44 day (p^ Cancri d) planet is shown with the phased HET 
data. Two cycles are plotted. The signatures of Cancri b, Cancri d, and p^ Cancri have 
been subtracted from the observations. 
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Fig. 3. — The RV orbit of the 14.66 (p^ Cancri b) day planet is shown with the phased HET 
data, Cancri c, Cancri d, and p^ Cancri e are subtracted from the observations. 
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Fig. 4. — The RV orbit of the 2.808 day planet (p^ Cancri e) is shown with the phased 
HET data. The Cancri b, Cancri d, and p^ Cancri d planets are subtracted from the 
observations. The rms of the HET residuals in this fit are 5.4 m s~^ . 



-16- 



30 p 1 1 1 1 1 1 — I — I — I 1 1 1 1 1 1 — I — 

25 - - 
20 - - 
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Period(days) 

Fig. 5. — Periodogram of the velocity residuals of the HET data, after the orbits of 
Cancri b, Cancri c, and p^ Cancri d have been subtracted from the original velocities. 
The peak at 2.808 days has a false alarm probability of 1.730796e-09. 




Fig. 6. — The Keplerian RV orbit of the 4 planets (combined) around p^ Cancri is shown 
with the HET data. 
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Table 1. Measured Velocities for Cancri 



Julian Date Radial Velocity Uncertainty 



m s ^ m s ^ 



o 


245292 / .yo4iUU 


48.3959 


A /1 

4.42 


.ug 20 


245292 ^9ooo9o 


45.2ii9 


4.oi 


2452920.9 i o95i 


32.9(^45 


A 00 

4.o3 


2452920.9 o555 


33. o595 


A nn 
4. (JO 


< 


2452929.9 /oiU i 


i0.i059 


A n 1 


m 


2452929.9 ( o099 


0.020i 




24529oi.9o4ioi 


-45.2540 


n^ 
3.96 


> 
in 


24529oi.9oD m4 


/in o /I /I n 
-49.244U 


A nT 
4.U 


24529o2.9DOoUi 


- /5.3o0o 


A no 
4.U3 


oo 
vn 


24529o2.9Doo94 


- / 4.000U 


/I 1 


00 

o 


24529o4.94oUd 


1 n 1 1 /I 
-iUi.i5o4 


A 

4.o5 


245z9o4.945 / io 


1 nn n/^n/^ 

-i(j(j.yo9o 


/I 01 


O 


2452945. 9iU r 22 






)-ph 


245294o.9ioo2 / 


'70 '71 rr 

-/o. loo 


A on 
4.0U 


o/itxon/in nonno/1 
2452949. 92UUo4 


or [r777 
-00.0 / / 


0.0 ( 




2452949. 922oo9 


QQ KA on 


A 00 

4.ZZ 




245295o.92i9Do 


5.5979 


4.93 


> 


2452953.924579 


5.7871 


4.59 




2452954.905716 


23.7456 


4.52 


:ar: 


2452954.908321 


21.0795 


4.10 


2452958.888656 


49.5613 


5.51 




2452958.893853 


56.6997 


6.23 




2452958.896458 


51.7819 


6.15 




2452959.910350 


7.8077 


4.47 




2452959.912966 


5.8177 


4.51 




2452963.873048 


-77.3188 


4.35 




2452963.875664 


-82.3257 


4.60 




2452965.859597 


-71.6258 


5.20 




2452965.862201 


-70.6408 


5.88 




2452981.832542 


-45.6285 


5.51 
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Table 1 — Continued 



Julian Date Radial Velocity Uncertainty 



2452981. 


.835147 


-43.6642 


5.22 


2452983. 


.819870 


6.0030 


5.70 


2452983. 


.822475 


6.4444 


5.14 


2452984. 


,824483 


35.3138 


5.76 


2452984. 


,827087 


35.9807 


5.40 


2452987. 


,025269 


70.3006 


5.51 


2452987. 


,027861 


58.3073 


5.35 


2452988. 


,026683 


46.7548 


5.16 


2452988. 


.818139 


28.2357 


5.88 


2452988. 


.820743 


27.0745 


5.61 


2452991. 


.026634 


-20.4009 


5.34 


2452991. 


,029250 


-20.5040 


5.38 


2452993. 


,794697 


-55.7202 


5.83 


2452993. 


,797301 


-61.1158 


5.95 


2452994. 


,789284 


-59.4158 


5.51 


2452994. 


,791889 


-61.9043 


5.19 


2452995. 


,021521 


-46.0055 


5.32 


2452995. 


,024125 


-45.6422 


5.33 


2452995. 


,801048 


-29.6404 


5.30 


2452995. 


,803641 


-29.0286 


5.00 


2452998. 


,021804 


29.4389 


6.51 


2452998. 


,784278 


45.7845 


7.14 


2452998. 


,789811 


42.9436 


6.40 


2452999. 


,012045 


41.6519 


5.15 


2452999. 


,014650 


41.9007 


5.09 


2452999. 


,769043 


56.6176 


5.55 


2452999. 


,771613 


51.6558 


5.54 


2452999. 


,997250 


57.7691 


5.46 


2453000. 


,775776 


70.4790 


5.09 


2453000. 


,778381 


71.2492 


5.53 
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2453003.787771 


35.7763 


5.84 


2453003.790375 


34.3674 


5.94 


2453003.999645 


34.7542 


5.59 


2453004.002261 


29.9420 


5.44 


2453006.766660 


-54.2333 


5.92 


2453006.769264 


-48.2144 


5.84 


2453024.931112 


-64.2770 


7.00 


2453024.933716 


-68.9715 


6.74 


2453031.911171 


69.6003 


7.00 


2453031.913775 


68.8923 


7.09 


2453033.683626 


13.3570 


6.78 


2453034.927281 


-11.3642 


6.37 


2453034.931622 


-10.2679 


6.65 


2453039.910567 


-30.3113 


6.99 


2453042.670770 


47.1661 


6.54 


2453042.673376 


50.6863 


6.76 


2453042.883142 


65.7295 


6.61 


2453042.885746 


66.7046 


6.59 


2453050.652384 


-41.9300 


6.95 


2453051.641174 


-47.8081 


7.29 


2453051.643778 


-53.3926 


7.86 


2453052.856436 


-69.5633 


6.17 


2453053.631560 


-57.1619 


6.84 


2453053.634164 


-59.2720 


6.83 


2453056.837680 


27.0449 


6.72 


2453057.620820 


49.7831 


6.04 


2453057.623424 


49.1894 


5.87 


2453063.596764 


-4.9079 


26.52 


2453063.839021 


-13.6337 


7.29 


2453063.841683 


-12.6510 


8.42 
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m 


m 


2453063.844287 


-16.8885 


17.02 


2453074.806906 


76.7781 


5.89 


2453074.809510 


76.9143 


5.71 


2453082.782947 


-48.3146 


5.24 


2453082.785551 


-46.4042 


5.00 


2453084.780423 


-11.6753 


4.92 


2453084.783027 


-10.9265 


5.28 


2453089.742959 


68.9187 


6.29 


2453089.758629 


67.9285 


5.87 


2453089.761221 


68.1677 


5.10 


2453095.747987 


-60.8753 


5.02 


2453095.750591 


-60.0757 


4.92 


2453098.760397 


-24.0996 


5.65 


2453098.762954 


-29.2928 


5.55 


2453101.724356 


50.2003 


4.71 


2453101.726960 


46.5767 


4.79 


2453102.719393 


66.5337 


4.55 


2453102.721997 


67.5490 


4.41 


2453103.714695 


66.9446 


4.61 


94^;'^ 1 (T? 71 79QQ 


fizL p;q'?q 


ZL ^P, 


2453110.698488 


-61.0881 


4.45 


2453110.701092 


-65.0179 


4.87 


2453111.689262 


-73.3320 


4.63 


2453111.691877 


-70.9320 


4.87 


2453112.687743 


-57.5881 


4.71 


2453112.690347 


-53.7139 


5.07 


2453116.667708 


61.5199 


6.12 


2453116.670312 


61.3680 


5.03 


2453117.697852 


62.8956 


5.02 
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